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Abstract: Aiming at the distributed flexible job shop scheduling problem under dual resource constraints con- 
sidering the influence of workpiece transportation time between factories and machines, a distributed flexi- 
ble job shop scheduling problem (DFJSP) model with the optimization goal of minimizing completion time 
is established, and an improved mayfly algorithm (IMA) is proposed to solve it. Firstly, the mayfly position 
vector is discrete mapped to make it applicable to the scheduling problem. Secondly, three-layer coding 
rules of process, worker and machine is adopted, in which the factory selection is reflected by machine 
number, according to the characteristics of the model, a hybrid initialization strategy is designed to improve 
the population quality and diversity; Thirdly, an active time window decoding strategy considering trans- 
portation time is designed for the worker-machine idle time window to improve the local optimization 
performance of the algorithm; In addition, The improved crossover and mutation operators is designed to 
expand the global search range of the algorithm. Finally, through simulation experiments, the results of 
various algorithms are compared to verify the effectiveness of the proposed algorithm for isomorphism and 
isomerism factories instances. 


Keywords: dual resource constrained; distributed flexible job shop scheduling; transportation time; 
improved mayfly algorithm; discrete mapping 


1. Introduction 


Economic globalization and the rapid development of new technologies such as 
cloud computing, big data and the Internet of things are leading the manufacturing in- 
dustry to the modes of production globalization, cloud manufacturing and intelligent 
manufacturing [1-4]. The traditional centralized manufacturing mode can no longer meet 
the fierce market competition and diversified customer needs. A large number of enter- 
prises expand their production to the distributed production environment. Through the 
rational allocation, optimal combination and sharing of resources of enterprises or facto- 
ries in different regions, they can quickly improve production efficiency and quality and 
reduce costs and risks [5-6]. Providing reasonable and efficient resource allocation and 
scheduling services for collaborative production of distributed factories is a research 
hotspot in the field of scheduling. 

Distributed flexible job shop scheduling problem (DFJSP) has many constraints and 
is challenging to solve. It is a more complex NP-Hard problem, which has attracted the 
attention of many scholars at home and abroad in recent years. De Giovanni et al.[7] pro- 
posed an improved genetic algorithm to solve DFJSP with multiple flexible manufactur- 
ing units, which extended the decoding of FJSP to include the information of FMU alloca- 
tion, and used a new local search operator to improve the population quality; Ziaee[8] 
assumes that the job shop of each factory/unit is configured as a flexible job shop, in order 
to obtain a high-quality scheduling plan quickly, a fast heuristic algorithm based on the 


constructive process was proposed for DFJSP; Sang et al.[9] built a DFJSP model for col- 
laborative manufacturing of smart factories, proposed a high-dimensional multi-objective 
memory algorithm combing with improved NSGA-III and local search methods, and op- 
timized economic indicators and green indicators; Xu et al.[10] proposed a three-layer 
coding of hybrid genetic tabu search algorithm to optimize the completion time, cost, 
quality and carbon emissions of DFJSP considering outsourcing some workpieces; Meng 
et al.[11] proposed four MILP models and a constraint programming (CP) model for 
DFJSP, and verified the effectiveness and superiority of the model through small and large 
instances. The above studies usually assume that a workpiece can only be processed in 
one factory. In recent years, some scholars have studied the open shop scheduling envi- 
ronment considering that a workpiece can be processed in different factories. Luo et al.[12] 
established a DFJSP mathematical model considering the workpiece transfer with the op- 
timization objectives of the completion time, factory load and energy consumption, and 
designed the GLS initialization method and a variety of neighborhood structures to im- 
prove the memetic algorithm for solving it; Gong et al.[13] considered that workpieces can 
be transferred between machines, job shops and factories, studied the distributed produc- 
tion scheduling of different factories and job shops, and proposed a new memetic algo- 
rithm to optimize the problem's completion time and total energy consumption; Du et 
al.[14] established a MILP model of dual-objective DFJSP considering the constraints of 
crane transportation and energy consumption, and proposed an EDA-VNS hybrid algo- 
rithm to solve it. However, there is still lack of relevant literatures on the transportation 
time required for workpiece transfer. 

In addition, most studies on DFJSP only consider the constraints of machinery and 
equipment, but ignore the constraint of worker resources that is inseparable from produc- 
tion scheduling. There are some related studies on FJSP in the integrated manufacturing 
environment. Meng et al.[15] studied dual resource-constrained flexible job shop sched- 
uling problem(DRCFJSP) with energy consumption awareness, proposed two MILP mod- 
els and designed a neighborhood search algorithm with eight neighborhood structures 
for the solution. Obimuyiwa[16] considered a limited number of cross-trained skilled in- 
stallation operators, established a detailed MILP model of DRCFJSP, and solved it with a 
genetic algorithm; Gong et al.[17] proposed a DRCFJSP model considering workers’ flex- 
ibility and proposed a hybrid artificial bee colony algorithm to solve it. 

Previously, our research group studied DRCFJSP of integrated manufacturing mode, 
which better solved the differences of the operation levels of workers, but ignored the 
influence of transportation time when in the model established [18]. Therefore, to further 
study the distributed manufacturing mode scheduling problem, this paper considers the 
influence of worker-machine dual resource constraints and workpiece transportation time, 
and constructs a distributed flexible job shop scheduling problem under dual resource 
constraints(DFJSPD) model with the completion time as the optimization objective. And 
an improved Mayfly Algorithm (IMA) is proposed according to the characteristic of the 
model, where discrete mapping of mayfly position variables, a new three-layer coding 
method for multiple resource constraints, a mixed initialization strategy, an active time 
window decoding algorithm for idle time window and the improved crossover and mu- 
tation modes are designed to improve solving performance of the algorithm. Finally, the 
experimental results show that the proposed DFJSPD model is more in line with the actual 
scheduling situation and the effectiveness and superiority of the improved algorithm. 

The remainder of this work is organized as follows. Section 2 establishes the DFJSP 
model and Section 3 gives the improved mayfly solution. Experimental analysis and re- 
sults are presented in Section 4. Section 5 concludes the main findings. 


2. DFJSPD Modeling 
2.1. DFJSPD Problem Description 


DFJSPD is described as: ” workpieces NSL) needs to be processed in F flex- 
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ible job shop type factories / . The workpiece ! has” processes with prior- 
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ity constraints among them. In the factory ’/, there are” machines al r) 


and wokas ” ("=L Zaat] operating the machines, generally’ 7 Wf, The ma- 
chines in the factory include CNC machines and non-CNC machines. The workers are 
responsible for the preparation including the loading and unloading, the replacement of 
tool fixtures and cleaning of CNC machine, and the preparation and machining of non- 
CNC machine. Individual differences among workers lead to the differences of operation 


efficiency. The basic times of preparation and machining on different machines are known. 


Workpieces are transported between machines in a factory or different factories with 
known transportation time. DFJSPD can be divided into four sub-problems: factory selec- 
tion, machine selection, worker selection and process sequencing. 

Taking the data in Table 1 as an example, the machines with asterisks are CNC ma- 
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AM Wa) Scheduling diagram of DFJSPD as shown in Figure 1, each color represents 


a job, which order by time course is the process number. 


Table 1. An example of DFJSPD. 
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Figure 1. Schematic diagram of DFJSPD scheduling. 


The relevant assumptions for scheduling are as follows: 

(1) All workpieces, machines and workers are available at time 0; 

(2) At any time, a machine or worker process one process at most; 

(3) Atany time, each workpiece is processed by one worker operating one machine 
at most; 

(4) Workers can only flow in the factory and the transfer time is considered within 
the preparation time and the loading and unloading time of the workpiece is 
considered within the transportation time; 
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When the front process on the machine where a process is located is the front 
process of the same workpiece, the preparation time is negligible. 

The transportation time of workpieces transferred between factories is longer 
than that transferred within factories; 

The transportation time of each workpiece before and after the first process is 
negligible; 

There are enough transport tools to complete the transfer of the workpiece; 

No interruption is considered in the processing. 


2.2. Mathematical modeling 


Based on the relevant descriptions and assumptions in Section 2.1, the DFJSPD sched- 
uling optimization model is constructed to minimize the maximum completion time 
(makespan). The parameters are shown in Table 2. 


Table 2. Parameter description 


parameter explanation parameter explanation 
i,h index of workpieces k,l index of machines 
I8 index of processes Sir index of workers 
k ; 
0O; j the front process of process O; on the machine k 
Ta basic machining time of process O; on machine k 
Tix basic preparation time of process O; on machine k 
Eks efficiency of worker $ operating machine k 
iiks actual processing time of process O; by worker S operating machine k 
S ijks starting time of process O; processed by worker S operating machine k 
Cits completion time of process O; processed by worker S operating machine k 
C; and C, completion times of factory f and workpiece i 
TY, Transportation time between machines k and | 
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X ijks 1, if process O; is processed by worker $ operating the machine k ; 0, otherwise 
x Lif i'=i and j'=j-1 „that is the front process the machine k where process O; is located is the 
X jij 
front process of the same workpiece ; 0, otherwise 
Yu 1, workpiece i is transported between machines k and | ; 0, otherwise 
H, 1, machine k is CNC-machines; 0, otherwise 


min(makespan) = min(max C, ) = min(max C,) 
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Equation (1) indicates that is the objective function is the maximum completion time 152 
of all factories, that is the maximum completion time of all workpieces; Equation (2) indi- 153 
cates that the actual operation time of workers is equal to the ratio of the standard time to 154 
the efficiency of workers operating the machine; Equation (3) indicates that the whole 155 
processing process is continuous without interruption, and the completion time of the 156 
process is the sum of the starting time, transportation time, actual preparation time and 157 
actual machining time; Equation (4) indicates that the actual start time of the process is 158 


restricted by process, transportation time, selected worker and machine factors; Equation 
(5) indicates that each process can only be processed on one machine by one worker at the 
same time; Equation (6) indicates that each workpiece can only be transported between 
two machines at the same time; Equation (7) indicates that the time of two processes pro- 
cessed on the same machine cannot be crossed; Equation (8) indicates that the time of two 
processes processed by the same worker cannot be crossed.. 


3. Improved mayfly algorithm for DFJSPD Model 
3.1. Design of Improved Mayfly Algorithm 


Mayfly algorithm(MA) is a swarm intelligence algorithm, which is inspired by the 
flight and mating behaviors of mayflies[19]. MA that combines the advantages of particle 
swarm optimization (PSO)[20], genetic algorithm (GA)[21] and firefly algorithm (FA)[22] 
is often used to solve continuous problems[23-24]. But DFJSPD is a discrete problem, so 
this paper improves MA to make it more suitable for the field of job shop scheduling. The 
improved strategies are as follows: 

(1) The positions of mayflies are discrete mapped to obtain the corresponding dis- 
crete codes of the feasible solutions, which make MA suitable for job shop sched- 
uling problem; 

(2) The mixed population initialization strategy, which designs a variety of initiali- 
zation methods based on the heuristic rules of time to improve the population 
quality and diversity; 

(3) Design active window decoding algorithm for the three-layer codes to obtain a 
better scheduling scheme; 

(4) Improve the crossover and mutation operators, which increase the population 
diversity and improve the global exploitation and local exploration of the algo- 


rithm. 
The flow chart of the improved mayfly algorithm is shown in Figure 2. 
Update velocities and positions 
of may flies 
Initialize parameters and velocities| | (1) Discretize the positions of mayflies 
and positions of mayflies | 1 to three-layer codes 
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: three-layer codes ; (4) Perform crossover and mutation 


ee ye ee TA operations of three-layer codes 
| (3) Decode and evaluate the target | eee ime ce eres 
| value of individuals | Merge offspring and parent 
: individuals. 
Update individual and 
global optimal value 


Whether the 
termination 
conditions are met? 


Evaluate all individuals and 
select male and female 
populations for next generation 


Output optimal results 


Figure 2. Flow chart of improved mayfly algorithm. 
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3.2. Discrete mapping of mayfly 


In the MA algorithm, the initialization formula of velocity and position of mayfly is 
as follows: 


P=(u,-1,)*rand() +l, (9) 


In formula (9), ” and l, represent the upper and lower boundaries of velocity or po- 


ta nx max(n, 
sition respectively. P is a set of random number vectors of (n) 


the initial velocity or position vector of the mayfly. 

The workpieces information cannot be read directly from the position vector, so it 
needs to be discrete mapped and converted into the process code of integer sequence. 
4 max(n,)=3 


, which represents 


Taking the example in Table 1, ” = , the mapping process is shown in Table 
3. In table 3, Pin is that Pm is arranged in ascending order. Ais the original index of Pr 
corresponding to Pm ; ¥2 is% /max(n,) rounded up, that is converted to workpiece number; 
gi is the occurrence order in turn of the corresponding workpiece in X , namely, the pro- 
cess number; So is that all invalid processes in gi are set to zeros; a set of workpiece num- 


bers in 22 corresponding to nonzero numbers in $2 is the process code OC. 


Table 3. Discrete mapping of the position vector of mayfly. 


index Pm P pa de Wa G So oc 


1 0.0402 
0.2580 7 3 1 1 3 
2 0.0613 
0.6934 11 4 1 1 4 
3 0.0977 
0.0977 3 1 1 1 1 
4 
0.6920 0.2580 1 1 2 2 1 
5 0.5983 
0.5983 5 2 1 1 2 
6 0.6301 
0.7139 8 3 2 2 3 
7 0.6919 
0.0402 10 4 2 2 4 
8 0.6920 
0.6301 4 2 2 2. 2 
9 0.6934 
0.7003 2 1 3 0 
10 0.7003 
0.6919 9 3 3 3 3 
11 0.7139 
0.0613 6 2 3 3 2 
12 0.9185 
0.9185 12 4 3 0 


3.3. Coding and Decoding 


According to the problem characteristics of DFJSPD, each mayfly individual repre- 
sents a solution that contains multiple information of process, factory, machine and 
worker. This paper adopts a three-layer coding method: process coding (OC), machine 
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coding (MC) and worker coding (WC). The three vector elements are one-to-one corre- 
spondence, and the factory selection reflected by machine coding. Taking the data in Table 
1 as an example, the coding is shown in Figure 3. The three-layer codes are read from left 


to right at the same time. Each number in OC represents the workpiece number, which J 


-th occurrence in turn and / is the process number of the workpiece. The corresponding 
positions of MC and WC represent the machine that processes the process and the worker 
who operates the machine. 


OC 


Optional 
machines 


Optional 1 
workers i 2 


Figure 3. Schematic diagram of coding. 


A feasible solution code can only show the resources allocation information, and the 
complete scheduling scheme needs to be obtained by decoding. The most commonly used 
decoding method is the insertion decoding, as shown in Figure. 5. Each process is ar- 
ranged in the earliest processing time of the available idle time window. This method can 
effectively reduce the waste of idle time, but the utilization of idle time window may not 
be high. 

Therefore, this paper designs an active time window decoding algorithm considering 
transportation time for DFJSPD. The key is that when a process can be properly adjusted 
within the available idle time window, the start time is determined at the latest end time, 
and the larger idle time window is left as far as possible for subsequent processes, as 
shown in Figure 5. According to the process, transportation time, machine type and flex- 
ible constraints of machine and worker in DFJSPD, whether the preparation time and ma- 
chining time of each process are affected is judged and arranged in an appropriate idle 


time window. The flow chart of decoding is shown in Figure 6, a aaa et wea ak 


indicates the start and end time of the idle time window. The main steps are as follows: 

Traverse three-layer codes to obtain a process, the selected machine and worker, and 
the start and stop time of the idle time window of machine and worker; 

According to formula (2)-(8), the actual processing time and the earliest start time of 
the process are calculated; 

Traverse the available idle time window, compare the remaining idle time window 


after insertion and “"¢7« to determine the actual start time of the process. ave, is the av- 
erage processing time of all processes on the machine, as shown in (10); 
ps X (X ngr Th gk + X nek Tek ) 
aveT, _ hal g=l — 
X > X hgk 
h=1 g=1 (10) 


Determine the actual end time of the process, update the idle time window, and de- 
code the next process until all the processes are completed. 
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Figure 5. Active time window decoding 


3.4. Mixed initialization 


Initializing the population is the premise of intelligent algorithm optimization, the 
initial population quality has an important influence on the convergence speed and opti- 
mization ability of the algorithm. When setting the initial solution, it is necessary to im- 
prove the quality of the solution according to the optimization objective. In addition, the 
diversity of the initial solution should be ensured to avoid the population falling into local 
optimum. In this paper, for OC that initial position vector of mayfly is discrete mapped, a 
variety of initialization rules considering time (population ratio is 4:2:2:1:1) are designed 
for the corresponding MC and OC. 
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(3) 


(4) 


Heuristic rule based on completion time. Each process in OC is traversed suc- 
cessively to select the worker-machine combination that can minimize the com- 
pletion time; 

Global initialization rule. Each process of each workpiece is traversed succes- 
sively, and the worker-machine combination for the shortest total cumulative 
processing time is selected for it. 

Local initialization rule. Each process of each workpiece is traversed in turn to 
select the worker-machine combination for the shortest cumulative processing 
time of the workpiece. 

Single initialization rule. Each process in OC is traversed successively to select 
the worker-machine combination that can minimize the processing time. 
Random initialization rule. Each process in OC is traversed successively to select 
any worker-machine combination that can process the process. 
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Figure 6. Flow chart of an active time window decoding algorithm 


3.5. Updating of mayfly 
Male mayfly in MA constantly adjusts its velocity and direction according to its his- 


torical optimal position (pbest) and the optimal position in the population (gbesi) to move 
towards the optimal position. When the male mayfly is in the best position, the wedding 
dance mode is executed, and its velocity varies according to inertia weight and dance co- 
efficient. The update formulas of the velocity and position of male mayfly are as follows: 


8 *myv; +a” ( pest, -x1,)+a,e"" ( ghest, =x; ) if f(x) > f (gbest) 
t . 
gy, tar, if f(x) < f(gbest) 


t+1 


(11) 


pee ee mv; (12) 


t . 


where, ™i is the velocity of the male mayfly/ in dimension / at time step’. £ is the 
inertia weight. ^ and ® are positive attraction constants used to scale the contribution of 


the cognitive and social component respectively. P is a fixed visibility coefficient. d is the 


nuptial dance coefficient. “is a random value in the range [-1,1]. "> is the Cartesian 
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formula (12), “iis the position of the male mayfly Í in dimension / at time step! . 

In order to reproduce offspring, female mayflies are attracted by male mayflies. They 
adjust their velocity and direction to fly towards male mayflies and keep approaching. 
when a female mayfly is not attracted by a male mayfly, so it flies randomly. The update 
formulas of the velocity and position of female mayfly are as follows: 


fit = ei +a,e°™ (xi my if FO) > FO) m 
g* fu, + A*r, if fy) fx) 


y =y t (14) 


Fi is the velocity of the female mayfly ! in dimension / 


where, at time step! . Pl is a ran- 


domly walk coefficient and ’nt is the Cartesian distance between male and female mayflies. 


In formula (12), Yi is the position of the female mayfly Í in dimension J at time step! . 

In order to balance the global exploitation and local exploration in the early and late 
stages of the algorithm, the dynamic inertia weight is adopted[19], which linearly de- 
creases with the number of iterations. After the mayfly position is updated, it is mapped 
to process code according to the discrete method in section 3.2. The original worker-ma- 
chine combination of each process is unchanged, and their positions changed accordingly. 
After updating the position of mayfly corresponding to the process code in Figure 3, the 
three-layer codes after discrete mapping are as shown in Figure 7. 


Figure 7. the codes after updating of mayfly. 


3.6. Crossover and mutation operators 


The crossover operation of the original MA algorithm is based on the quality of the 
solution. A male mayfly breeds with and the female mayfly at the same fitness level with 
it. But the original crossover formula is applicable to the continuous optimization problem 
[23]. Therefore, in this paper, crossover and mutation operators are designed respectively 
for the characteristics of the three-layer codes. 

(1) Crossover operator 

(a) The IPOX crossover operator is adopted on OC. (b) An improved IMPX crossover 
operator [24] is adopted on MC. If the machine is unavailable after exchanging, a machine 
with the shortest processing time is selected from the set of alternative machines for the 
corresponding process. (c) Since it is assumed that workers cannot transfer across factories, 
a two-point crossover operator considering factory constraint is designed for WC. For the 
parent P1, two crossover points are selected and exchanged respectively with any two 
workers in the same factory in the parent P2. If the worker after the exchange is unavail- 
able, a worker with the highest efficiency is selected in the optional workers set for the 
corresponding machine. Three crossover operators are shown in Figure 8, where P repre- 
sents the parent and C represents the offspring. 
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Figure 8. Crossover operators. 


(2) Mutation operator 

(a) Reverse sequence mutation is adopted on OC. The genes reverse sequence ex- 
change between two different mutation sites randomly selected and MC and WC are ad- 
justed accordingly to keep the original worker-machine combination of each process un- 
changed. (b) The mutation rules on MC is the same as OC, and if the machine after mu- 
tating is unavailable, the worker-machine combination is reselected in the optional ma- 
chines and workers set. (c) Two points of mutation is adopted on WC that is selecting 
randomly two different workers in a factory to exchange positions. If the worker after the 
exchange is unavailable, a worker with the highest efficiency is selected in the optional 
workers set for the corresponding machine. The three mutation operators are shown in 


Figure 9. 
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(b) Mutation of MC (c) Mutation of WC 
Figure 9. Mutation operators. 


4. Experiment and Analysis 
4.1. Testing instances and parameter setting 


Since there are few previous studies on DFJSPD, there is no the relevant benchmark 
instances for reference and test. This paper design 20 DFJSPD instances of isomorphism 
and isomerism factories by extending 10 FJSP benchmarks of Brandimarte[27] that are 
widely applicable in the field of FJSP and used by a large number of researchers [9, 28, 29]. 

(1) Isomorphism factories instances (SMk01-SMk10) 

The quantity and functional flexibility of machines and workers in isomorphism fac- 
tories are the same. Set the workpieces processing time and the worker-machine of infor- 
mation and copy these to construct factories F1 and F2 with the same processing environ- 
ment. 

(2) Isomerism factories instances (DMk01-DMk10) 

There are differences in the number and function flexibility of machines and workers 
in isomerism factories. Split the machines information and the number of workers of fac- 
tory F1 in isomorphism factories instances to construct factories F1 and F2 respectively, 
and workers’ skills of each factory are complete flexibility. 


For any instance, the process preparation time is generated in CRA the worker 
efficiency is generated randomly in Ja , the transportation time between machines 


U[1,5] 


in the same factory is generated in , and the transportation time between machines 


in different factories is generated in U[6,10] [9]. 30% of the machines are selected as CNC 
machines. The scale and machines information of each instance are shown in Table 4. 


nx (mx w) x F represents ” workpieces are processed in F factories, and” machines and 
W workers are in each factory. 


Table 4. Descriptions of the DFJSPD instances. 
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Isomorphism Scale CNC Isomerism Scale 
factories factories CNC machines 
instances nx (mxw)x f machines instances nx(mxw)x f 
SMk01 10x(6x4) x2 2,4,8,10 DMk01 10x(3x2) x2 2,4 
SMk02 10x(6x4) x2 2,4,8,10 DMk02 10x(3x2) x2 2,4 
SMk03 15x(8x6) x2 1,8,9,16 DMk03 15x(4x3) x2 1,8 
SMk04 15x(8x6) x2 1,3,9,11 DMk04 15x(4x3) x2 1,3 
SMk05 15x(4x3) x2 3,7 DMk05 15x(2x2) x2 3 
SMk06 10x(15x8) x2 2,3,4,7,17,18,19,22 DMk06 10x((8x4) + (7x4)) 2,3,4,7 
SMk07 20x(5x4) x2 2,7 DMk07 20x((3x2) + (2x2)) 2 
SMk08 20x(10x6) x2 1,3,10,11,13,20 DMk08 20x(5x3) x2 1,3,10 
SMk09 20x(10x6) x2 2,4,8,12,14,18 DMk09 20x(5x3) x2 2,4,8 
SMk10 20x(15x8) x2 2,3,7,10,17,18,22,25 DMk10 20x((8x4) + (7x4)) 2,3,7,10 
To verify the superior performance of IMA, comparing with the results of traditional 
MA, GA, PSO and FA algorithms. Traditional algorithms adopt three-layer encoding, ran- 
dom initialization and insertion decoding methods. The parameters setting of the algo- 
rithms are shown in Table 5. All experiments are run on MATLAB R2020a. 
Table 5. Parameters setting of algorithms 
Number of 
Algorithm Population scale Other Parameters 
Iterations 
female:50, a =1,a, =a, =15,f=2,d= fl =1, 
IMA 200 
male:50 Emax = 1-3, 8 min = OF 
female:50, 
MA ave 200 a, =1,a, =a, =1.5,8 =2,d = fl =1,g =0.8 
male:50 
GA 100 200 Crossover rate 0.8, mutation rate 0.05 
PSO 100 200 Inertia weight @ = 0.8, learning factor €, =C, = 2 
attractiveness of firefly 4) =1, light absorption 
FA 100 200 


coefficient Y =!, random parameter & = 0.3 


4.2. Results analysis 


In order to avoid random differences, each instance is solved 20 times by each algo- 
rithm, and the results of isomorphism and isomerism factories instances are shown in 


min ave 


Figure 6-7. max and Cinax are the optimal and mean value of 20 results of an algorithm. 
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Table 6. The results of isomorphism factories instances 


IMA MA GA PSO FA 

Cis C imax C C max Cien C max Cii C imax Cna C max 
SMk01 42.82 43.86 45.21 46.26 46.83 47.66 46.49 47.38 47.22 48.04 
SMk02 31.78 32.39 32.93 33.38 33.10 33.43 33.70. 34.66 34.53 35.11 
SMk03 154.28 156.06 155.38 157.77 155.86 156.89 164.85 165.28 163.59 165.03 
SMk04 63.83 65.63 67.60 68.39 66.51 67.89 68.53 69.99 67.70 69.51 
SMk05 142.02 145.27 147.78 149.69 149.21 156.52 148.93 154.08 150.30 158.20 
SMk06 91.66 93.02 93.94 94.77 94.12 95.09 96.18 97.57 95.69 96.47 
SMk07 122.00 124.35 126.06 128.36 126.30 128.19 128.67 131.26 127.30 130.64 
SMk08 430.00 435.96 445.59 447.21 441.00 447.72 446.53 453.16 447.31 451.48 
SMk09 319.41 321.98 329.40 332.58 330.24 334.33 349.96 351.68 342.18 346.66 
SMk10 237.13 240.55 245.76 256.69 248.64 259.62 252.78 256.63 260.37 262.41 

Table 7. The results of isomerism factories instances 
IMA MA GA PSO FA 
ew ge a a e. O a re 
max max max max max max max 
DMk01 68.03 69.97 70.03 73.27 74.28 76.63 74.33 76.08 75.79 77.00 
DMk02 45.56 47.38 46.03 47.99 46.54 48.63 48.74 50.60 48.39 50.70 
DMk03 235.26 257.21 268.34 284.60 263.26 279.46 266.92 287.95 276.92 286.13 
122.00 
DMk04 104.51 111.52 114.19 117.74 112.84 118.54 118.45 : 116.93 122.60 
Cmax 

DMk05 234.53 243.02 256.25 270.02 256.93 268.72 258.51 271.35 260.77 272.13 
DMk06 133.89 142.71 153.49 158.90 150.54 159.85 153.72 161.22 157.08 160.73 
DMk07 213.13 225.70 228.69 238.58 232.32 239.04 249.30 253.57 249.30 252.71 
DMk08 735.70 765.51 788.57 811.98 770.23 805.96 801.65 836.58 781.34 844.57 
DMk0o9 529.02 547.75 569.38 585.13 550.38 572.43 556.30 585.96 571.10 591.50 
DMk10 428.64 447.54 458.24 463.94 446.04 469.25 468.68 480.76 466.86 479.06 


Table 6-7 shows that due to the randomness of population initialization and evolu- 
tion process, there are differences among the results of each example. Overall, the optimal 
value and mean value of IMA are less than other algorithms. The larger the scale of the 
example, the more obvious the difference. It shows that the active time window decoding 
can better coordinate resources and reduce idle time, especially when solving large-scale 
problems. 

To significantly compare algorithms differences, two evaluation indicators are de- 
signed according to literature [29]: the relative percentage deviation of the optimal value 
MRPD and the relative percentage deviation of the mean value ARPD. The calculation 
formulas are as follows: 
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MRPD = x100 


low 
max 


(15) 


s i _ plow 
ARPD = Lys, 100] 
S j=1 


(16) 


In formula (15) ~ (16), for any instance, Cmax is the ! -th result of an algorithm solv- 
low 


ingit, “max is the optimal value of all the results of all the algorithms, 5 is the number of 
solving. the values of MRPD and ARPD are shown in Table 8-9. 


Table 8. Evaluation results of isomorphism factories instances. 


badness te 7 w IMA MA GA PSO FA 
MRPD ARPD MRPD ARPD MRPD ARPD MRPD ARPD MRPD ARPD 

SMk01 42.82 0 2.44 5.58 8.04 9.36 11.31 8.57 10.65 10.28 12.19 

SMk02 31.78 0 1.91 3.62 5.04 4.15 5.20 6.04 9.06 8.65 10.47 

SMk03 154.28 0 1.16 0.71 2.26 1.02 1.69 6.85 7.13 6.03 6.97 

SMk04 63.83 0 2:77. 5.86 7.10 4.15 6.31 7.31 9.59 6.01 8.84 

SMk05 142.02 0 2.29 4.06 5.40 5.06 10.21 4.87 8.49 5.83 11.39 

SMk06 91.66 0 1.49 2.49 3.39 2.68 3.74 4.93 6.45 4.39 5.24 

SMk07 122.00 0 1.93 3.33 5.21 3.52 5.07 5.47 7.59 4.34 7.09 

SMk08 430.00 0 1.39 3.63 4.00 2.56 4.12 3.84 5.39 4.03 5.00 

SMk09 319.41 0 0.80 3.13 4.12 3.39 4.67 9.56 10.10 7.13 8.53 

SMk10 237.13 0 1.44 3.64 8.25 4.85 9.48 6.60 8.22 9.80 10.66 

Table 9. Evaluation results of isomerism factories instances. 
ETET s v w IMA MA GA PSO FA 
MRPD ARPD MRPD ARPD MRPD ARPD MRPD ARPD MRPD ARPD 

DMk01 68.03 0 2.85 2.94 7.70 9.19 12.64 9.26 11.84 11.41 13.18 
DMk02 45.56 0 4.00 1.03 5.34 2.15 6.73 6.98 11.07 6.21 11.28 
DMk03 235.26 0 9.33 14.06 20.97 11.90 18.79 13.46 22.39 17.71 21.62 
DMk04 104.51 0 6.71 9.26 12.66 7.97 13.43 13.34 16.73 11.88 17.31 
DMk05 234.53 0 3.62 9.26 15.13 11.19 16.03 10.22 15..70 9.55 14.58 
DMk06 133.89 0 6.59 15.39 18.83 12.44 19.39 14.81 20.41 17.32 20.04 
DMk07 213.13 0 5.90 7.30 11.94 9.00 12.16 16.97 18.98 16.97 18.57 
DMk08 735.70 0 4.05 7.19 10.37 4.69 9.55 8.96 13.71 6.20 14.80 
DMk09 529.02 0 3.54 7.63 10.61 4.04 8.20 5.16 10.76 7.95 11.81 
DMk10 428.64 0 4.41 6.91 8.24 4.06 9.47 9.38 12.16 8.92 11.76 


From Table 8-9, it can be seen that all the MRPD values of IMA for 20 DFJSPD in- 
stances are 0, which indicates that the optimal results of IMA are the best in all algorithms. 
In most instances, the differences of MRPD values between MA and GA are small, and 
the probability of obtaining the suboptimal solution is greater than PSO and FA, which 
indicates that the crossover and mutation operations of MA and GA can expand the search 
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range and improve the probability of obtaining the better solution, while the operation of 
updating the position according to the better solution of PSO and FA has low solution 
accuracy and is easy to fall into local optimum. Comparing with ARPD values, the average 
performance of IMA is better than MA, GA, PSO and FA, which shows that IMA hybrid 
initialization strategy and improved crossover and mutation operators improve the solu- 
tion quality and robustness of the algorithm, as shown in Fig. 10. 


90 


85 


oo 
© 


Makespan 


x 
a 
` 
z 
z 
= 
è 
L 


70 


0 50 100 150 200 
Generation 


Figure 10. Algorithms iteration process diagram of DMk01. 


The convergence performance of the algorithm is analyzed by taking example 
DMk0O1. It can be seen from Figure 10 that the hybrid initialization strategy of IMA im- 
proves the quality of the initial population. The five algorithms can converge to their re- 
spective optimal solutions, but the convergence speed of GA and FA is slow and the so- 
lution accuracy is low converging to 74.28 and 75.79. PSO has fast convergence speed, but 
it is easy to premature converging to 74.33. The convergence speed of IMA and MA is fast 
in the early stage, and the completion time converges to 70 or so in the 30-th generation. 
The wedding dance and random flight operations in the late iteration make the algorithm 
jump out of local optimum, but the accuracy of the crossover and mutation methods of 
MA is not high converging to 70.03. The multiple crossover and mutation operators of 
IMA can expand the search range, find optimal solution 68.03. Overall, the optimization 
ability and speed of IMA are better than other four algorithms. 

From the Gantt charts of the optimal solutions of SMk01 and DMkO1 in Figs. 11-14, it 
can be seen that the scheduling schemes of the same number of workpieces in distributed 
factories with different structures are different, and the final completion time is also dif- 
ferent. Workpieces are circulated among multiple open structure factories, and the pro- 
cessing resources of each factory can be fully utilized in time. The workpieces are trans- 
ferred among multiple open structure factories, which can make full use of the processing 
resources of each factory timely. The more comprehensive resources, the higher the pro- 
cessing efficiency. This paper can provide more clear and accurate scheduling scheme for 
distributed factories collaborative intelligent manufacturing with different structures. 
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Figure 11. The Machine Gantt chart of SMkK01 with Optimal solution. 426 
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Figure 12. The Worker Gantt chart of SMk01 with Optimal solution. 428 
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Figure 13. The Machine Gantt chart of DMk01 with Optimal solution. 430 
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Figure 14. The Worker Gantt chart of DMk01 with Optimal solution. 


5. Conclusions 


In the actual distributed flexible job shop production process, the worker-machine 
dual resource constraints and the transportation time of workpiece flowing among ma- 
chines seriously affect the scheduling plan. In order to construct a scheduling model that 
is more in line with the actual processing situation, this paper establishes the DFJSPD 
model with the optimization objective of the minimum of maximum completion time. For 
the multiple constraints of the discrete DFJSPD problem, Mayfly algorithm is improved 
by using the discrete mapping method of continuous variables, three-layer coding, decod- 
ing algorithm of an active time window, hybrid initialization strategy and the crossover 
and mutation operator of three-layer codes to improve the global exploitation and local 
exploration performance of MA. Finally, the basic instances are expanded to design the 
instances of isomorphism and isomerism factories. The simulation experiments are car- 
ried out and compared with various algorithms. The results show that IMA's solving qual- 
ity and speed are significantly better than those of other algorithms. The model and algo- 
rithm in this paper can help enterprises realize distributed collaborative intelligent man- 
ufacturing. 

However, the above research is still insufficient. It needs to be further optimized, 
such as only optimizing single objective, without considering the impact of transport re- 
source constraints, and without considering resource changes in a dynamic real-time 
scheduling environment. Therefore, this paper will also study the dynamic scheduling of 
DFJSP under multi-resource constraints, optimize economic indicators and green indica- 
tors such as energy consumption and noise, and explore intelligent algorithms with better 
performance in the future. 
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